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j-». , A rigorous method for introducing the variational principle describing relativistic ideal hydrody- 

namic flows with all possible types of breaks (including shocks) is presented in the framework of an 
exact Clebsch type representation of the four-velocity field as a bilinear combination of the scalar 
Ch , fields. The boundary conditions for these fields on the breaks are found. We also discuss the local 

3 ■ invariants caused by the symmetries of the problem, including relabeling symmetry. In particular, 

the generalization of the well-known nonrelativistic Ertel invariant is presented. 
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i^, a. Introduction. In the paper we discuss some problems related to the ideal relativistic hydrodynamic (RHD) 

1-^ ' flows in the framework of the special relativity. They are pertinent to description of the flows with breaks in terms 
I , of the canonical (Hamiltonian) variables based upon the corresponding variational principle, and introducing the 
^ ' generalization of the Ertel invariant. These subjects are of interest both from the general point of view and are 
'^-H I very useful in solving the nonlinear problems, speciflcally, for the nonlinear stability investigation, description of the 
c/2 . turbulent flows, etc. The necessity to consider the relativistic flows is motivated by a wide area of applications to the 
^ ^ ■ cosmological problems. 
C/3 . Variational principles for the ideal relativistic hydrodynamic (RHD) flows are widely discussed in the literature, 

^^' see, for instance, [l|, y, y| and citations therein. As for the nonrelativistic flows, the least action principle is conve- 
nient to formulate in terms of the subsidiary fields and corresponding velocity representation known as the Clebsch 
representation, see 0, 0, S S S S Il3l ■ These subsidiary fields can be introduced explicitly by means of the Weber 
transformation, [ll|, see also 0, I3- Alternatively, they naturally arise from the least action principle as Lagrange 
multipliers for necessary constraints. Using these variables allows one to describe the dynamics in terms of canoni- 
cal (Hamiltonian) variables. The nontrivial character of the Hamiltonian approach is due to the fact that the fiuid 
dynamics corresponds to the degenerated case, see [13, UM ■ 

In the papers [l4L Il5l | it was shown that the hydrodynamic flows with breaks (including shocks) can be described 
1^^ ^ in terms of such least action principle, which includes (as natural boundary conditions) the boundary conditions for 
f^ ■ the subsidiary fields. In the nonrelativistic case the triplet of the subsidiary fields corresponds to the Lagrange labels 
■^ I of the fiuid particles, say, fj.^ , which are advected by the fiuid, 



dtH^^O, dt = dt + vV, B = 1,2,3, (1) 



where v denotes three-velocity. These equations along with the entropy advection and the fluid mass conservation 
are assumed as constraints. Corresponding Lagrange multipliers, As, 9 and (p, along with /i^ enter the Clebsch type 
^ '' velocity representation, 

^' /9v = -pVip - AbVm^ - evs, (2) 



1 J ■ where p and s denote fluid density and the entropy per unit mass. 
rS I b. Variational principle. The relativistic least action principle can be formulated in a close analogy to the non- 

jrt ■ relativistic one. Namely, introduce action A, 



f d^xL, 



A^ / d^xC, (3) 



with the Lagrangian density 



C^-e{n,S) + GJ"Q,a,, 

G=(l,z/B,e), Q^iip,p'',S), 5 = 1,2,3, 

where vb, ©, f, M^ represent subsidiary fields; n, S and e(n, 5) denote the particle's number, entropy and energy 
proper densities, J" = nit" is the particle current, and u" is the four-velocity, u" = u^{l,v/c), vP = l/i/l — v^/c^; 



comma denotes partial derivatives. Small Greek indexes run from to 3, and the Latin indexes run from 1 to 
3; x^ = ct, r = (x^,x^,x^). The metric tensor, g"'^, corresponds to the flat space-time in Cartesian coordinates, 
g°'0 = diag{ — 1, 1,1,1}. The four- velocity obeys normalization condition 

u"uq = gai3U°'u'^ = -1- (5) 

Below we consider the four-velocity and the particle density n as dependent variables expressed in terms of the 
particles current J", 



w" = J"/|J|, n^\J\ = v/-J"Ja. (6) 

The fluid energy obeys the second thermodynamic law 

de = nTdS + •nT wdn = nTdS + Wdn , (7) 

where T is the temperature and w = e +p is the proper enthalpy density, p is the fluid pressure, W ~ w/n. 

The action in Eq. Q depends (for a flxed or inflnite volume) on the independent variables J", 6, and Q = {ip, /i^, S), 
A = A[J", ip, fi^ , S', j^B , 6]. Its variation results in the following set of equations 

dV" :=^ Wuc, = Va^ -GQ^a , (8) 

S^ :=^ J% = 0, (9) 

Sfj,^ :=^ do,{J"vB) = 0, or Dua = 0, (10) 

5vB ■■=> Dn^ = 0, (11) 

5S :=> dc,{rQ), or DQ = -T, (12) 

dQ:=^ DS = 0, (13) 

where D = u^da- Eq. ||HJ) gives us Clebsch type velocity representation, cf. Ref. Q- Contracting it with u" results 
in the dynamic equation for the scalar potential ip, 

Dip = W. (14) 

Both triplets ji^ and vb represent the advectcd subsidiary fields and do not enter the internal energy. Therefore, it 
is natural to treat one of them, say, n^ as the flow line labels. 

Taking into account that the entropy and particles conservation are incorporated into the set of variational equations, 
it is easy to make sure that the equations of motion for the subsidiary variables along with the velocity representation 

sponds t 
fluid stress-energy- momentum T"^ conservation, cf. Ref. p^Tlll^. 

Qrpap 



reproduces the relativistic Euler equation. The latter corresponds to the orthogonal to the flow lines projection of the 

B 



dxP ~ ^^ 

T°'f = wu'^u'^ + pg'^'^ . (16) 

Note that the relativistic Euler equation can be written as 

(14,/j - T//3,a)w^ = TS^o. , (17) 

where the thermodynamic relation 

dp = ndW - riTdS (18) 

is taken into account. Vector Va, sometimes called Taub current, |la |. plays an important role in relativistic fluid 
dynamics, especially in the description of circulation and vorticity. Note that W can be interpreted as an injection 
energy (or chemical potential), cf., for instance [l3|, i.e., the energy per particle required to inject a small amount 
of fluid into a fluid sample, keeping the sample volume and the entropy per particle 5* constant. Therefore, Va is 
identified with the four-momentum per particle of a small amount of fluid to be injected in a larger sample of fluid 
without changing the total fluid volume and the entropy per particle. 



c. Boundary conditions. In order to complete the variational approach for the flows with breaks, it is necessary 
to formulate the boundary conditions for the subsidiary variables which do not imply any restrictions on the physically 
possible breaks (the shocks, tangential and contact breaks), are consistent with the corresponding dynamic equations 
and thus are equivalent to the conventional boundary conditions, i.e., to continuity of the particle and energy- 
momentum fluxes intersecting the break surface R{x°') = 0, cf. Ref. [l8l |. 

{J} = 0, J = J"na , (19) 



{T'^%} = 0, (20) 

where Ua denotes the unit normal vector to the break surface. 



na = NcjN, N^^Rc. N=^/nJ^, (21) 

and braces denote jump, {X} = X\ji=^q — X\]^-^o. 

Our aim is to obtain boundary conditions as natural boundary conditions for the variational principle. In the 
process of deriving the volume equations we have applied integration by parts to the term J"GdQ_a. Vanishing of the 
corresponding surface term along with that resulting from the variation of the surface itself lead to the appropriate 
boundary conditions after the variational principle has been specified. 

Rewriting the (volume) action with the break surface being taken into account in the explicit form as 

A= f d*xY^ 'C'^^(?i?), (22) 

where 9 stands for the step-function, we obtain the residual part of the (volume) action in the form 

^^\res ^ I d^^Y. ['i^^D{R)SR + e{<;R)do,{rG5Q)] . (23) 

Here 5d denotes Dirac's delta-function and we omit index ? labeling the quantities that correspond to the fluid regions 
divided by the interface i? = 0; superscript ? ^ corresponds to the quantities in the regions i? ^ 0, respectively. 
Integrating the second term by parts and supposing that the surface integral J d'^x X]r=±i ^a {0{(;R){u"GSQ)) vanishes 
due to vanishing of the variations 5Q at infinity, we arrive at the residual action expressed by the surface integral 

^Kes = f <^^^Y. '^^DiR) [CSR - R,a.J"GSQ\ . (24) 

•' ';=±l 

SQ here means the limit values of the volume variations, 5Q = {5Q)fi-±o. It is convenient to express these variations 
in terms of variations of the boundary restrictions of the volume variables, 6{Xfi-±o) = 6X , and variation of the 
break surface. It is easy to show that 

'SX = 6X+ \N\-'^n°'X^a6R ~ X^^P^f^Sf'^ , (25) 

where P"^ — (5"^ — n^-np. and 5f^ is an arbitrary infinitesimal four- vector related to the one-to-one mapping of the 
surfaces i? = and R + 5R = 0. 

Vanishing of the action variation with respect to variations of the surface variables 5R and 5f^ (which are supposed 
to be independent) results in the following boundary conditions 

5R -.^ {p + {u'^n^fw} = 0, (26) 



Sff^ :^ P'^p {WTNaU^} =^ 0, or P'^^ [JWu^] = 0, (27) 

which are equivalent to continuity of the momentum and energy fluxes, cf. Eq. 120|l . Here we consider that the 'on 
shell' value of the volume Lagrangian density, Ceq, is equal to the pressure, Ceq = — e + nGDQ = —e + w ~ p. 

Now we can complete formulation of the variational principle appropriate both for continuous and discontinuous 
flows. The independent volume variables are indicated above, and independent variations of the surface variables are 



6R, Sf^, variations of the surface restrictions of the gencrahzed coordinates Sip, Sfi^ , supposed to be equal from both 
sides of the break, {6(p} = {Sn^} = 0, and 6S with {6S} ^ 0. Under these assumptions we arrive at the following 
subset of the boundary conditions 

8^ :^ { J"n„} = { J} = for {5^} = 0, (28) 

S'jl^ :^ {vBrria} = J{vb] = for {SJi^} = 0, (29) 

bS^ -^ rnj^^ = je± = 0. (30) 

Eqs. (|26|I - H28() reproduce the usual boundary conditions, and Eqs. H29|l . (|30|l arc the boundary conditions for the 
subsidiary variables. Other boundary conditions for the latter variables do not strictly follow from the variational 
principle under discussion. But we can find them from the corresponding volume equations of motion, providing, for 
instance, that they are as continuous as possible. |22| The natural choice corresponds to continuity of their fluxes, 

{nau°'rni^} = J{ii^} = 0, (31) 

{nau°'nip] = J{lp} = 0. (32) 

The set of the boundary conditions given by Eqs. (PH|) " (j^ is complete and allows one to describe any type of breaks, 
including shocks. For the latter case J ^ and we arrive at continuity of the variables vb, A*^, <P and zero boundary 
value of Q. For J = the flow lines do not intersect the break surface and we obtain very weak restrictions on the 
boundary values of the subsidiary variables, cf. nonrelativistic case discussed in Refs. [ij^ [l3- Note that for the 
specific case J = (slide and contact discontinuities) we can simplify the variational principle assuming all both-side 
variations of the subsidiary variables to be independent. 

The above variational principle allows modifications. First, it is possible to exclude constraints, expressing the 
four- velocity by means of representation ||SJ|. In this case the volume Lagrangian density can be chosen coinciding 
with the fluid pressure, cf. Ref. |3|, where the continuous flows arc discussed in detail. Second, we can include into 
the action the surface term respective for the surface constraints, cf. Refs. |l4. flSl Il9l l20| . where such surface terms 
are discussed for ideal hydrodynamics and magnetohydrodynamics in the nonrelativistic limit. This can be done for 
the cases both with excluded and non excluded volume constraints. 

Canonical variables. Starting with the action of Eq. Q and Lagrangian density given by Eq. Q) we can 
introduce the canonical (Hamiltonian) variables according to the general receipt. Namely, let Q represents the 
canonical coordinates. Then 

SA 

P=j-- = .fG = i7r^,n^s,ns) (33) 

gives us conjugate momenta. Relations ()33f) cannot be solved for t he g eneralized velocities Qq suggesting that we are 
deahng with the degenerated (constraint) system, cf. Refs. |3, la, Il3i 113 ■ -^ut the constraints are of the first type. 
Thus, performing the Legendre transform with respect to Q we arrive at the Hamiltonian density 

H = PQ,o-p{W,S), (34) 

where we suppose that the four- velocity is given by representation JHJ. Making use of the definition (|33() and of the 
time component of the velocity representation, Eq. ^ , we can transform the first term in Eq. H34(l as 

PQfi = J°GQ,o = -TT^^^o = TTvV°. (35) 

Taking into account the normalization condition for the Taub current, VaV°' = —W^, we obtain 



1/° = y/W^ + VaV" . (36) 

Consequently, we arrive at the following Hamiltonian density 

n = H{P, Q, 0,a; W) = VW^ + VaV' t:^ - p{W, S). (37) 



In terms of the canonical coordinates and momenta the space components of the velocity are 

TT^K = -PQ.a- (38) 

The canonical equations following from this Hamiltonian reproduce in a 3 + 1 form the above dynamic equations for 
the variables entering the Taub current representation. Variation of the action with respect to the chemical potential 
W results in 

(39) 



v/1 + VaV^/W^ 

Obviously, this relation is equivalent to Eq. (|36|l , expressing the particle density n in terms of the variables entering 
the Hamiltonian. 

Underline, that the Hamiltonian given by Eq. (|37|1 depends not only on the generalized coordinates (p, fj.^ , S, their 
spatial derivatives and conjugate momenta, but also on the chemical potential W as well. Evidently, we can consider 
W as the additional generalized coordinate with zero conjugate momentum, tt\y = 0. This condition is consistent 
with the dynamic equations due to the fact that Oottw = dTi/dW = 0, cf. Eq. (|39|l . 

Bearing in mind the flows with breaks one can see that in the above discussed variant of the least action principle 
we do not arrive at the additional surface variables except that defining a break surface, R. But it enters the 
action functional without derivatives. Therefore, corresponding conjugate momentum is zero-valued. Introducing the 
Hamiltonian variables for the flows with breaks we have to treat R as the surface function, defining some (surface) 
constraint. The latter is nothing else than continuity of the normal component of the fluid momentum flux, Eq. ]2(i\i . 

d. Poisson brackets. The Poisson brackets in the canonical variables are of a standard form. Symbolically, 

{Q-^ix),P[,iy)}^5iS{x~y), (40) 

where S{x — y) is spacetime Dirac's delta, x and y are spacetime points, Q is shorthand notation for (p, fi and S, 
and, analogously, Pg denotes corresponding conjugate momenta, tt^, tt^b and irg. 

e. Ertel invariant. In addition to energy, momentum, and angular momentum conservation, for the ideal hy- 
drodynamic flows there are exist specific local conservation laws related to the dragged and frozen-in fields, and 
corresponding topological invariants (vorticity, helicity, Ertel invariant, etc.), cf. Refs. 0j S S 113 123 ^'^'^ citations 
therein for the nonrelativistic case. They are caused by the relabeling symmetry. Discussion of these problems for 
the relativistic flows seems insuflficient, see Refs. [llSSIlfl and citations therein. Exploitation of the above descrip- 
tion permits one considering these invariants to be simplifled. For example, consider here generalization of the Ertel 
invariant for the relativistic fluids (to my best knowledge, this item was not discussed earlier). Defining the Ertel 
four-current, 

1 

^PuS^y = -*uj"''S^^ , (41) 

Here e"'^'*'^ is Levi-Civita tensor, w/j^^ is the (Khalatnikov) vorticity 
Moreover, the Ertel four-vector £" is proportional to the particle 

£"=EJ°', (42) 

in view of £"^ = resulting in E = £°/ J° being dragged by the fluid, 

DE = 0, (43) 

i.e. E is the scalar invariant of the motion. In the nonrelativistic limit it coincides with the Ertel invariant (curlv • 
Vs)/p, where p denotes the fluid density. 

/. Helicity current. The helicity invariant in the nonrelativistic case exists for the barotropic flows and presents 
pseudoscalar v • curlv. The strict analog for the relativistic case is the pseudovector 

Z^ = If^^'^LOp^y, = *Lo''''V, , (44) 

Strict calculations show that for the isentropic flows the helicity current Z" is conserved, Z"" ^ = 0. 

For the general type flows there exists generalization of the helicity current. Namely, consider reduced Taub vector, 

Va,=Va, + e^,„ , (45) 
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where obeys Eq. (|12() . and the corresponding reduced vorticity tensor, 

i^ap = Vf3,a - Vq,,/3 . (46) 

This tensor is orthogonal to the flow lines, 

^apu^ = 0, (47) 

and the reduced helicity current 



is conserved for arbitrary flows. 



Z" = *S"''K (48) 



,Q = W K,Q = ^ ^ ^ai. = 0. (49) 



g. Conclusion. We have shown that it is possible to describe the relativistic ideal fluids with all physically 
allowable breaks in terms of the least action principle both in the Lagrangian and Hamiltonian description. The 
boundary conditions for the subsidiary variables entering the Clebsch type velocity representation are obtained in 
two different ways: one part follows from the variational principle as natural boundary conditions while the other one 
was obtained from the dynamic equations under assumption relating to absence of the corresponding sources and the 
maximal continuity compatible with the volume equations. Note that it is possible to change the variational principle 
in such a way that all boundary conditions will result from it, i.e., they become natural boundary conditions. For this 
purpose it is necessary to modify the variational principle by addin g a surface term with corresponding constraints, 
similarly to the nonrelativistic case (compare with the papers [3 uM for the hydrodynamics and [l^ |23| for the 
magnetohydrodynamics). This variants are to be discussed in the forthcoming papers. 

The approach discussed allowed us to give a simple treatment of the additional invariants of the motion, in particular, 
to present generalization of the Ertel invariant for the relativistic flows. This approach is suitable for the general 
relativity and for the relativistic magnetohydrodynamics as well. Note that for the flows without breaks the general 
relativity case is discussed in detail in the paper |2] ■ The discontinuous flows for the general relativity can be described 
in analogy to the above discussion and will be published elsewhere. 
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